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Abstract 

We considered the Inverse Compton Scattering (ICS) of charged par- 
ticles onto photons whose distribution is a Black Body Radiation (BBR) 
deriving the exact energy and angular differential distribution in the gen- 
eral case and in its most useful expansions. These results can be success- 
fully applied in high energy accelerators experiments to evaluate the ICS 
contribution from the thermal photons in the cavity as well as in astro- 
physics where the ICS of cosmic rays plays a relevant role in a variety of 
phenomena. In particular we show how our formulae reproduce the ICS 
energy spectrum recently measured at LEP, how it could be considered 
a key tool in explaining the Gamma Ray Bursts (GRB) , SGRs energy 
spectrum. Finally we predicted the presence of a low gamma flux, nearly 
detectable at hundred of TeV from SNRs SN1006 as well as, at lower en- 
ergy (tens TeV, due to gamma ray cascading in cosmic BBR), from relic 
extragalactic highest cosmic rays sources born by jets in AGN,as blazars 
3C279,Mrk421,Mrk 501. 
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Introduction 



The Inverse Compton Scattering (ICS) plays a relevant role in highest en- 
ergy astrophysics (cosmic rays and gamma astronomy) [1-6] and high energy 
physics (LEP I, LEP II, accelerators) [7-9]. Indeed from one hand the ICS of 
high relativistic cosmic rays (either electron at GeV or proton and nuclei at 
much higher energies) onto electromagnetic fields (either cosmological Black 
Body Radiation (BBR) at T 2.73 K, interstellar lights, radio waves or even 
stationary magnetic fields) is the source of high energy photons (X, gamma 
rays) which we do observe in the Universe as diffuse or point source, on the 
other hand the ICS is often the main process responsible for the slowing down 
(i.e. for the energy losses) of energetic charged particles; indeed ICS is often the 
main cause of cosmic rays lifetime behaviour as well as of their energy spec- 
trum depletion at harder regions (i.e. of their detailed spectrum shape and 
evolution) [5-6]. Moreover, ICS of relativistic jets by compact objects onto 
thermal photons by a star companion in a binary system or accretion disk 
might be, as recently proposed [10], the key process able to produce "gamma 
jets" responsible (by their rotation and blazing into different directions) of the 
puzzling Gamma Ray Bursts (GRB). So even if this subject seems apparently 
settled [3-4] we revisit the ICS process in order to obtain an analytic and com- 
pact formula able to describe the differential energy and angular ICS scattered 
photon number spectrum of relativistic charges onto BBR photon spectrum. 
One of the main feature of our final expressions is that we can easily cover the 
entire range of energy of the ICS energy spectrum. These results improve the 
Montecarlo simulation because the latter must consider only a thin portion 
of the energy range in order to inspect in detail the ICS spectrum [9]. In the 
following we describe our approach and we show some remarkable results. 

The ICS onto BBR spectrum 

We follow a standard procedure to get the general ICS spectrum of a relativis- 
tic charge (electron, proton, nuclei in cosmic rays or in accelerators bunches) 
hitting photons whose distribution is a BBR spectrum. We consider first the 
photon target distribution in the Laboratory Frame (LF) where the BBR is 
isotropic and homogeneous, then we transform it to the Electron Frame (EF) 
where the BBR is still homogeneous but highly anisotropic, therefore we eval- 
uate in the EF the usual Compton scattering and finally we transform back 
the diffused differential photon number to the LF. The starting photon target 
distribution, in the LF, is the well-known isotropic and homogeneous BBR 
whose number density per unit energy e Q and solid angle Q Q is given by the 
Planck formula 
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where Kb is the Boltzmann constant and T the BBR temperature. We trans- 
form this distribution to the EF by standard Lorentz boosts choosing as z 
axis the direction coincident with the initial electron momentum and remind- 
ing that dn /de is a relativistic invariant [2]. In the following we label by a 
* the quantities related to the electron frame EF, by a subscript if they are 
considered before the scattering, by a subscript 1 if after; so we have 

cos9* +(3 # 
° = 1 + (3cos9* ' (p o = <Po, e = 7e (l + (3cos0 o ) 

where f3 is the adimensional electron velocity and 7 the corresponding Lorentz 
factor. The transformed BBR number density distribution exhibits, in this rel- 
ativistic limit and in the EF, a clear dipole anisotropy and it is becomes more 
and more peaked around 9* = n as 7 increases. An analogous dipole signa- 
ture, at non relativistic regime, is the cosmological one found at millikelvin 
level in the T = 2.73 K BBR due to the Earth motion. The associated energy 
spectrum is mostly "green shifted" showing a maximum around 7/t_eT with 
respect to the original "red" BBR spectrum with a maximum around k^I '. 
The next step is to derive the total number of diffused photons in the EF; this 
number can be obtained as follows: 



dNf dn* Q da c 



dt\dt\dVt\de* dVt* de*dtt* de\dVt\ 
where d ^^* is the Compton differential cross section 
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c is the speed of light and the scattering angle must be expressed as a function 
of the other angles involved, i.e. the incoming ip* and the outcoming 9*, </?* 
angles 

cos#* c = sin 9* sin 9\{comp* cos ip\ + sin ip* a sin ip\) + cos 9* Q cos 9\ . (4) 

The last step is to derive the final exact ICS differential number distribution 
in the Laboratory Frame (where the scattered photons will be observed) by 
using the inverse Lorentz transformations; the result is: 
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where from eq.4 cos6> sc must be expressed as a function of Q*,0\,<p* o ,<p\. This 
"blue shifted" distribution is different from the previous ones because now the 
expected peak around a given e\ energy has been spread into a wide plateau 
from kbT up to 7 2 «bT energies. The most general ICS differential distribution 
can be obtained by means of a numerical integration of eq.5 over Q* but as 
we are more interested in high energy phenomena we also show how to sim- 
plificate the above formula in this cases. First of all we remind that in particle 
accelerators we are dealing with ultrarelativistic particles, i.e. 7 > 1; more- 
over for most astrophysical problems the photon energy in the EF is much 
smaller than the electron rest mass, i.e. e* <C mc 2 ; so it is often possible to 
approximate the Compton differential cross section by the Thomson one; as 
a third further approximation we may consider the ultrarelativistic-Thomson 
limit where the two previous conditions 7 ^> 1 and e* <C mc 2 are both satis- 
fied. Let us discuss these three different expansions. In the first case (7 ^> 1) 
the BBR photons, in the EF, are pratically all incident head-on so the incident 
angle 9* Q can be approximately written as 9* Q m n — K Consequently the scat- 
tering angle is related to the 9{ angle by the simple formula 9* c + 9\ ~ n and 
cos#* c ~ — cos 6^ within I/7. Moreover the kinematics of the ICS shows that 
the scattered radiation, in the LF, is strongly concentrated in a narrow cone 
Oi ~ 1/7 along the direction of motion of the ultrarelativistic particle. The 
ICS differential distribution in eq.5, in the ultrarelativistic expansion, reduces 
to an analytical expression 
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where 1 — f3 ~ I/27 2 and d£li : 
is the Thomson limit e* <C mc 2 
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: 2n9id9 1 . The second possible approximation 
In this case we can neglect all terms of order 



e*/mc with respect to 1 in eq.5 and the resulting differential distribution 
becomes 
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Indeed, from the previous eq.6-7, we can get the last expansion, the Thomson- 
ultrarelativistic formula. It stems from assuming 7 S> 1 and e* <C mc 2 at the 
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same time. This means that the ICS differential distribution can be obtained 
by setting e\ <C mc 2 7 in eq.6 and by setting 7 3> I in eq.7. In the first 
derivation we obtain: 



dN x 2-KK B Tr 2 c , 
-ei In 



dhdeM c^h 3 ^ 2 
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In the second derivation the integral contained in eq.7 can be simplified re- 
membering that 0* o ~ 7r — ^ so sin 9* cos 9* ~ I/7 and the integrand is sig- 
nificatively different from zero only inside a thin cone whose aperture is just 
I/7. Thus the result is a formula differing from eq.8 only by a factor of or- 
der 4j and in our assumptions this term is completely negligible. So the still 
analytical eq.8 is our ultrarelativistic-Thomson distribution. In fig. 1 we show 
the surprising metamorphosis of a Planckian Black Body into a taller (by a 
suppression factor 7~ 2 ) serial of smooth truncated hills at larger relativistic 
regimes. This process is mainly due to the overlapping of a series of blueshifted 
spectra obtained at different angular directions; the softer photons are the 
more isotropic ones while the harder photons are strongly anisotropic and lo- 
cated in the inner cone of the beam. The " Rayleigh" regions of the differential 
"BBR" spectra overlap each other even if calculated at different 9\ angles 
(0 < 9\ < it) while their peaks are higher and more and more blueshifted for 
9\ angles approaching zero. Consequently the exponential decay of the ICS 
spectrum on the right side reflects the "Wien" regions behaviour of the angu- 
lar "BBR" spectra at 9\ « 0. We give here a qualitative summary of the ICS 
energy spectrum behaviour in the whole e\ energy range. For e x <C K^T/7 2 
the spectrum exhibits a linear growth (like the Rayleigh region of a BBR 
one); for energies kbT/j 2 < e± < kbT the linear behaviour is modified by 
a logaritmic correction and the growth is proportional to e\ \n.{i^ 2 KbT / e\); 
for KbT < €1 < A^ 2 kbT the spectrum is spread into a very slowly linearly 
decreasing plateau up to e\ ~ 4^y 2 KbT; for e\ > A^ 2 k b T the spectrum decays 
as ei exp(— ei/A^kbT) (we remember the reader that the Wien region of a 
BBR distribution decays as e\ exp(— €i/k b T)). We remind that the commonly 
and widely applied analytical formulae on ICS are the ones derived by F.Jones 
[3] and based on ICS onto monochromatic and isotropic radiation; this ficti- 
cious and artificial "BBR" leads to a final spectrum in disagreement with the 
experimentally observed ones. In the next two sections we show that we can 
evaluate in the right way ICS spectra in high energy physics and astrophysics 
by means of our formulae. 



The ICS spectrum in LEP experiments 



We could directly verify the validity of present results by a (successful) fit 
of experimental ICS spectra obtained at LEP by A.Melissinos group [7] and 
by G.Diambrini-Palazzi group [8] (at a higher degree of precision). Indeed 
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the LEP vacuum pipe can be considered as a black body cavity at room 
temperature (T « 291 K); hence the electromagnetic radiation in thermal 
equilibrium is scattered and beamed ahead by the electron (and positron) 
bunches whose energy is E e = 45.6 GeV. At this energy corresponds a Lorentz 
factor 7 = 8.92 • 10 4 so we can apply our ultrarelativistic formulae for ICS 
to evaluate the interesting region of the spectrum. We compare our results 
with the Montecarlo simulations performed at LEP and able to fit the ICS 
effect in the experimental data. We compute either the Thomson and the 
Compton spectra in order to show the differences between the two curves and 
with respect to the Montecarlo simulation. The spectra can be obtained by 
a numerical integration over the 9\ angle in expressions 6-8 respectively for 
Compton or Thomson limit: 



where N e ~ 1.37 • 10 11 is the number of particles in a bunch and Ar = - ~ 
2 • 10~ 6 s is the flight time in the LEP straight section / ~ 600 m [9]. The three 
ICS spectra, derived from our equations and from the Montecarlo simulation 
are shown in fig. 2 and labelled respectively as: 

1) M: Montecarlo simulation (ref. [9] fig.6) 

2) T: Thomson approximation, from eq.8 

3) C: Compton approximation, from eq.6 

Note that our Compton spectrum shows a good agreement with the Monte- 
carlo simulation while the Thomson one, at high energies, is a factor 3 higher. 
This overestimate of the Thomson spectrum is clearly related to the inde- 
pendence of the Thomson cross section with respect to the photon energy. 
We notice that our ICS spectrum, in the Compton limit, is nearly coinci- 
dent with the Montecarlo approximation at low energy but for higher energies 
(> GeV) it is a 26% above. This discrepancy seems not to be related to our 
approximations because from our data we obtain, for the total event number, 
N\ = 2.646 a value pratically coincident (within 0.1%) with N-y = 2.65 found 
by Di Domenico [9] by Montecarlo simulations. The difference might be due 
to the Montecarlo method whose statistical procedure implies a smaller num- 
ber of events at higher energies. However this discrepancy does not affect the 
beam lifetime because in its evaluation only the number rate dNi/dti is in- 
volved. We remind that the ICS may also be used to study the bunch internal 
structure and the result of a coherent emission by the charges, in this case the 
optimal experimental set up is reached when the relativistic bunches are hit 
by collinear back photon emitted by a laser [11]. 



dNy(TH,KN) 

de\ 




f dNi(TH,KN) 

fix dtydeydQy 



dQyAr 



(9) 
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The ICS in Astrophysics 

The X — 7 astronomy traces mostly the presence of relativistic electrons 
(or, at lower level, of relativistic nuclei) by their synchrotron radiation or 
their ICS onto infrared, interstellar or cosmic radiation. Moreover the same 
ICS may become the main slowing down process for relativistic charges once 
magnetic field energy densities ps are below the corresponding cosmic photon 
energy densities Pbbr- This situation generally occurs in extragalactic spaces 
where ps <C Pbbr- In this framework cosmic rays electrons around SN1006 
have been recently discovered indirectly by their non thermal X-ray emission 
due to synchrotron radiation from ultrarelativistic electrons (7 > 10 s ) [12]. 
It is therefore of great interest and actuality to provide not just an order 
of magnitude for the corresponding ICS 7 ray flux but also for its detailed 
spectrum for such relativistic and ultrarelativistic cosmic rays electrons. We 
show in fig.3 these ICS spectra for 7 = 10 4 , 10 6 , 10 8 . It is important to note 
the existence of a high energy transition in the ICS from the Thomson to the 
Compton behaviour. This change occurs dramatically at the highest energy 
range of the spectra. The 7 Lorentz factor for Compton behaviour occurs for 
electrons (e~ + 7 — > e~ + e + + e~) and for protons (7 + p — >■ n + 7r + ,p + n°, 
...) in ICS with T = 2.73 K BBR respectively at huge energy values 

2m e c 2 m e c 2 9 2m 7T c 2 m^c 2 n 
7e = - - = — — ^ 1.6 - 10 y ; 7 P ~- — -=2.4-10 n 10 

< e Q > k b T < e G > k b T 

The ICS behaviour in the Compton limit may be qualitatively predicted keep- 
ing in mind that the ICS energy spectrum falls off at photon energies near 
€1 «3 7 2 /t^T and the energy conservation calls for a cut-off of the extension of 
the ICS spectrum plateau from kbT to mc 2 7 instead of reaching the usual ex- 
treme value 7 2 «bT. Let us better understand this ICS behaviour change from 
Thomson to Compton regime as follows: as long as the incident photon, in the 
EF, has a characteristic energy e* w ^k b T < mc 2 the diffused photon, in the 
EF, mantains most of its original energy and it is spread around nearly isotrop- 
ically. Therefore, once the spectrum is reviewed in the LF after the Lorentz 
boost, the previous different angular distribution of the photons in the EF be- 
comes, in the LF, a different energy distribution of the same photons. The an- 
gular integral of these differential "BBR" spectra dt ^ dVLl , whose "Rayleigh" 
regions overlap, produces in the energy range KbT < e\ < ^ 2 kbT the wider 
plateau shown in fig. 4. However, at Compton regime, when •ynsT > mc 2 the 
differential cross section a ~ ot ■ m a * \ leads, in the EF, to a diffused 
and anisotropic photon distribution which becomes more and more beamed 
in a "Compton cone" at angles Q* c < \j I fr^f\p^ ( n °t to be confused with the 
thinner "Lorentz" angles 9* L < ^ related to the boost transformations). All the 
photons diffused outside the Q* c cone can reach in the LF final energies of order 
ei w 7 2 /t B T(l-/3 cos ^) (1-/3 cos # s * c ) ~ 27 2 /t B T(l-/3cos^ c ) < mc 2 7. On the 
other side, around and inside the Compton cone Q* c but far from the edge of the 
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inner Lorentz cone 6* L (for example we may consider \\j ™ bT < @* sc < \j ™ b t 
where the majority (> 75%) of the scattered photons is found) the final 
avalaible energy, in the LF, will be e 1 ~ ^ 2 k b T(\ - (3cos6* sc ) > k b T + m §-'y. 
In conclusion the Compton behaviour of the ICS spectra, derived analitically 
from the exact formula 5 and described in fig. 5, piles up the photons near the 
edge energies in the range e\ < mc 2 7 and leads to an unexpected peak higher 

2 

and higher at those highest energy values as the Lorentz factor 7 ^> in- 
creases. It is also important to notice that this extreme Compton regime may 
arise in a different way: for example in non relativistic or relativistic cases 
when k b T ^> mc 2 . 



The main scenarios for ICS applications 



For a synthetic but complete picture of all the above described ICS behaviours 
we suggest to define the following characteristic regimes each labelled by the 
relevant parameters involved. 

1) (7 ~ 1, k b T <C m ^-) The non-relativistic Compton scattering limit 
onto cold- warm BBR (fig. 6). Here the simplest ICS spectrum becomes a 
self similar Planck spectrum; its reflectivity efficiency depends on the usual 
quantities n 7 , a T , ... One of the mostly celebrated application is the Sunyaev- 
Zeldovich effect in cosmology. 

2) (7 > 1, k b T < 2^) The ICS at relativistic Thomson limit onto 
cold- warm BBR (fig. 1-2). This is the case with most application in high en- 
ergy physics and astrophysics (LEP I, LEP II, GRBs, cosmic rays at energies 
E e < 10 14 eV and E p < 10 17 eV,...). Now the ICS spectrum deviates from 
a pure Planckian spectrum leading to a "cut-hill" spectrum with its smooth 
edge covering the energies from e± ~ k b T up to photon energies e± ~ 7 2 /teT. 
This " simple" ICS spectrum may be ruling the spectra of ICS by charge (elec- 
tron) jets (7 ~ 10 3 — 10 4 ) onto nearby stellar companion photons (e G ~ 0.5 eV) 
leading to successful gamma jets able to explain the integral spectra of GRBs 
[13]. This arguments have widely been developed by the authors and are still 
under consideration [14]. Because of the "smooth" behaviour of the edge of 
the spectrum there is not a one to one relation between the cosmic rays elec- 
tron or proton spectral index and the final ICS 7 rays spectral index. Another 
scenario where ICS leads to 7 rays takes place in extragalactic blazars which 
eject beamed cosmic rays (protons, nuclei,...) for huge distances. Their ICS 
onto 2.73 K BBR may also lead to collinear 7 jets; their energy is dominated 
by the electron presence (over the nuclei one) in the cosmic rays jets. Such 
gamma jet is the large scale version of the minijet model plus ICS considered 
by the authors for GRB production in galactic binary systems. 

3) (7 > 1, **f- < k b T < mc 2 ) The ICS at relativistic Compton limit 
onto cold- warm BBR (fig. 7). These energy windows are relevant for the 
highest cosmic rays interactions onto BBR [E e ^> 10 14 eV, E p > 10 20 eV). The 
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ICS spectrum exhibits a pile up of photons at the edge of the highest e\ energies 
leading to a peaked maximum at e\ ~ mc 2 7 and a sharp cut off at energies just 
above it. The presence of such a peak has important consequences in the cos- 
mic rays-7 rays links. The primordial incident cosmic rays spectrum at highest 
energies (power law, ...) leaves its original inprint into a similar "photocopy" 
7 ray spectrum. Moreover the new born high energies 7 rays (ei > 100 TeV) 
may also successfully interact with the BBR photons (electron pairs produc- 
tion) leading to electron-proton cascades or electromagnetic showers at lower 
and lower energies. The showers will arrest their growth as soon as the last 
degraded photon energies will reach the threshold E c = ^^mc 2 = 10 15 eV. 
It is important to remind that the 7 — 7 — > e + + e~ process at ultrarelativi- 
stic limit has a very similar cross section ( Breit- Wheeler, 1934) as the e + e~ 
annihilation (Dirac, 1930) or the Klein- Nishina (1929) cross section. One of 
the consequences is the expected cosmic background presence of relic TeV 7 
rays components (due to the direct cosmological ICS or due to its secondary 
showers at energies e\ < 100 TeV). This argument will be discussed in detail 
elsewhere. 

4) (7 ~ 1, k b T > mc 2 ) The non relativistic Compton scattering onto 
"hot" BBR (fig. 8). The resulting ICS spectrum deviates from the original 
one and exhibits a small peak at energies e\ ~ The effect as above is 
related to the energy dependence of the Klein-Nishina cross section. 

5) (7 > 1, mc 2 < k b T < mc 2 i) The relativistic ICS onto "hot" BBR. 
The spectrum shows a new "plateau" extending from e\ ~ k b T energies up 
to ei ~ mc 2 7 energies with a marked peak at e\ ~ mc 2 7 energies and a sharp 
cut off above, this behaviour is similar of point 3). 

6) (7 > 1, k b T = mc 2r )) The equilibrium regime of ICS onto "hot" 

BBR (fig.9). The spectrum shows a marked peak at e± ~ mc 2 7 energies with 
a sharp cut-off above. The known astrophysical scenarios where such ICS may 
play a role could be the earliest hottest (thermal equilibrium) cosmological 
epochs (t < 1 s, T > MeV) and in the hot thermal cores of supernovae 
explosions (k b T > 5 MeV). It is interesting to notice the nature of the non 
equilibrium ICS spectrum and it may be worthfull to derive the exact kinetic 
equations (due to such ICS) for the multicomponent thermal bath of the early 
universe as well as the SN explosive processes. Finally if the fireball model 
is a real event as the one needed to explain the GRB puzzle then such ICS 
(onto the last layers of the fireball explosion) would be smeared out by mul- 
tiscattering during last stages of the fireball into a final nearby thermal GRB 
spectrum. We do not recognize such a presence of thermal inprint in GRB 
[13]. 

7) (mc 2 < mc 2 7 < k b T) The ultrarelativistic "ultrahot" ICS (fig. 10). 
This ICS spectrum exhibits a peaked spectrum as usual at e± ~ mc 2 7 energies 
but also a significative decaying component at higher energies. We can see a 
"shoulder" which may become, in the extreme cases k b T 3> mc 2 7, a slowly 
decaying "plateau". Such a peculiar process (at the present) seems very hypo- 
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thetical but we have shown it in order to cover the entire range of possible ICS 
behaviours. Possible exotic scenarios where such extreme conditions (5,6,7) 
may take place efficiently occur also near miniblackhole evaporations in the 
early Universe where the cosmological BBR has a temperature lower or compa- 
rable or even greater than the corresponding earliest miniblackhole tempera- 
ture {nsTmbh ~ mc 2r y). This hybrid thermal bath, depending on miniblackhole 
primordial masses and distribution, may be dominant in non standard early 
cosmology bariogenesis and even in later cosmological nucleosynthesis. 



Conclusions and Applications 



The developed ICS formulae have a wide range of applications; in particular 
in fitting LEP I, LEP II experimental data and in understanding recent GRB 
puzzling spectra (in this last case we modified eq.8 in order to take into account 
the diluted and anisotropic BBR spectrum seen from the jet and we assumed 
a ring-like photon source ([14] eq.20). Moreover the extreme spectra of ICS at 
relativistic Compton limit onto cold-warm BBR, its peak at highest energies, 
may be probed at LEP I, LEP II using a diffused thermal optical light (a 
flash) in the beam pipe during the bunch crossing. Finally the recent evidence 
for cosmic rays electron at energies above hundreds TeV by their observed 
synchrotron radiation at soft X spectra implies the coexistence of a low but 
nearly detectable component of high 7 cosmic rays at energies E 1 < mc 2/ y w 
100 TeV. Because of arguments in 3) the cosmic rays electrons spectrum 
(dN/dE e ~ E~ 2 ) will be reflected also in the "photocopy" 7 spectrum. Their 
total energy flux 7 will be of the same order of magnitude of the X rays one: 

^7f^(^r <»> 

The corresponding 7 ray flux number will be extremely suppressed because of 
the X — 7 energy ratio 

dN 7 ^ 7 < E x > o in _ n/ r B E e 



dN x (j) X <E 1 > ^6fiGJ v 100 TeV 

and it follows that the needed area for such a low flux is possibly below the 
present air shower arrays sensibility. However an accurate direction-flux cor- 
relation might be able to observe in a near future a tiny 100 TeV flux. We 
suggest to seriously consider a detailed observational program in order to 
detect above the tiny 100 TeV 7 flux from SN1006 and possibly the other sug- 
gested candidates (Cas A, IC443, Tycho SN, ...) from SNRs. Their presence 
at the above fluxes is a necessary and compelling consequence of fundamental 
QED and cosmological BBR theories combined in present ICS models. Fi- 
nally we interpret the two X lobes around SN1006 as generated by an electron 
jet scattering onto the relic giant shell contrary to the more popular idea of 
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a Fermi shock acceleration mechanism. In these beaming models one should 
also expect a rare and strongly time dependent TeV 7 rays "burst" which 
could be observed in a short period of time (hours) once the observer is inside 
the thin jet cone direction, with an amplified integral intensity by a factor 
7 Lorentz larger then a diffused spherical source. Well known candidates are 
blazars or quasars as 3C279 , AGNs, as NGC 3079, or most recent TeV Mrk 
sources. Similar arguments leaded us to expect a low variable relic background 
of gamma (ten TeV) noise due to the pile up of cosmic integral ICS gamma 
rays and their electromagnetic cascades just below the electron pairs creation 
threshold on cosmic BBR. 
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Figure Captions 

Fig.l: The ICS spectrum from eq.8 for T=291 K and non relativistic 
(BBR) 7 = 1 (dot), and relativistic 7 = 10 (dot dash), 7 = 10 2 (dash), 
7 = 10 3 (continuous) Lorentz factor 

Fig. 2: The ICS spectra at LEP: Montecarlo simulation M (ref.9 fig. 6), 
Thomson approximation T (eq.8), Compton approximation C (eq.6) 
Fig.3: The ICS spectrum for T=2.73 K and 7 = 10 4 (dot dash), 7 = 10 6 
(dot), 7 = 10 8 (continuous) Lorentz factor 

Fig.4: The ICS "Thomson" spectrum for T=2.73 K and 7 = 10 4 Lorentz 
factor 

Fig. 5: The ICS spectrum evolution from "hill-like" Thomson spectra 
= 10~ 2 (dash), 10" 1 (lower continuous), 1 (dot dash) to Comp- 
ton peaked spectra = 10 (dot), 10 2 (higher continuous) for 7 = 10 
Lorentz factor 

Fig.6: The ICS Thomson spectra for T=291 K and 7 = 1 (BBR) (dot), 
7 = 3 (dash) Lorentz factor 

Fig.7: The ICS Compton spectra for k b T = 2.3 • 10~ 2 and 7 = 10 9 
Lorentz factor 

Fig. 18: The ICS Compton spectrum for k b T = mc 2 and 7 = 5 Lorentz 
factor 

Fig.9: The ICS Compton spectrum for k b T = 5.11 • 10 6 eV and 7 = 10 
Lorentz factor (k b T = mc 2r )) 

Fig. 10: The Compton ICS spectrum for k b T = 10 3 mc 2 (continuous), 
k b T = 10 4 mc 2 (dot) and 7 = 10 Lorentz factor 
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